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Self-similar networks with scale-free degree distribution have recently attracted much attention, 
since these apparently incompatible properties were reconciled in £jj by an appropriate box-counting 
method that enters the measurement of the fractal dimension. We study two genetic regulatory 
networks (Saccharomyces cerevisiae Q and Escherichai coli Q) and show their self-similar and 
scale- free features, in extension to the datasets studied by Moreover, by a number of numerical 
results we support the conjecture that self-similar scale-free networks are not assortative. From our 
simulations so far these networks seem to be disassortative instead. We also find that the qualitative 
feature of disassortativity is scale-invariant under renormalization, but it appears as an intrinsic 
feature of the renormalization prescription, as even assortative networks become disassortative after 
a sufficient number of renormalization steps. 

PACS numbers: 05.45.Df, 89.75.Hc, 87.16.Yc 



I. INTRODUCTION 

Until very recently, the celebrated properties of a scale- 
free degree distribution seemed to be incompatible with 
self-similar features of networks, in which the number of 
boxes Nb of linear size £b scales with £b according to 
a power-law Nb ~ £^ B ' with an exponent that is given 
by the fractal dimension ds- Using the box-counting 
method, Song et al. showed that many scale-free (SF) 
networks observed in nature can have a fractal structure 
as well Q. This result is striking, because the tiling and 
renormalization according to the linear size of the boxes, 
in which all pairs of nodes inside a box have mutual dis- 
tance less than Ib , appear to be physically relevant rather 
than being a formal procedure. Therefore, the essential 
quantity in the tiling is the linear size of the box, £b, 
defined by £b — 1 being the maximal distance between 
the nodes of the box. In the renormalization procedure 
the boxes are contracted to the nodes of the renormal- 
ized network whose edges are the interconnecting edges 
between the boxes on the original network. 

In this paper we study the genetic regulatory net- 
work of two well-known organisms, Saccharomyces cere- 
visiae and Escherichia coli We first determine 
the degree-distribution P(k), that is the probability for 
finding a node with degree k, to read off the exponent 7 
according to P(k) ~ fc -7 Q in order to check that the 
networks are scale-free. Next we measure Nb /N for vari- 
ous box-sizes £b to obtain the fractal dimension oIb from 
Nb /N ~ £' B dB . After renormalizing the networks accord- 
ing to the procedure proposed in [J , we measure the scal- 
ing behavior of the degree k' according to k' = s(£s)k, 
where k 1 stands for the degree of a node in the renor- 
malized network, k is the largest degree inside the box 
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that was contracted to one node with degree k' in the 
renormalization process, and s(£b) is assumed to scale 
like s{Ib) ~ £~B dk with a new exponent dk- The invari- 
ance of 7 under renormalization and the transformation 
behavior of the degree itself imply the relation p] 

7=1 + d B /d k (1) 

between the exponents. Therefore we check this relation 
by measuring ds , dk and comparing the values of 7 from 
Eq^with the measured 7 from the degree distribution. 

One of the important features of networks is their "de- 
gree" of assortativity. The notion of assortative mixing 
was known from epidemiology and ecology when it was 
introduced as a characteristic feature of generic networks 
by Newman (5J. Assortativity refers to correlations be- 
tween properties of adjacent nodes. One particular prop- 
erty is the (in- or out-)degree of a node as the number 
of its (in- or out-) going links, respectively. Degree-degree 
correlations can be recorded as histograms; in order to fa- 
cilitate the comparison between networks of different size, 
they can be also characterized by the Pearson coefficient. 
The Pearson coefficient is obtained from the connected 
degree-degree correlation function (jk) — (j) (k) after nor- 
malizing by its maximal value, which is achieved on a 
perfectly assortative network. Here, (jk) stands for the 
average of having vertex degrees j and k at the end of an 
arbitrary edge. The Pearson coefficient r takes values be- 
tween — 1 < r < 1, it is positive for assortative networks 
(r = 1 for complete assortativity) and negative for dis- 
assortative ones. We have measured this coefficient for a 
number of self-similar scale- free networks and present the 
results below. The reason why this feature is of interest 
in the present context is its relation to the power-law or 
exponential behavior of Nb(£b)- In particular, we are in- 
terested in the question whether disassortativity is scale- 
invariant on a qualitative level under renormalization ac- 
cording to the prescription proposed in pj , and why these 
properties go along. Disassortative features in protein in- 
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teraction networks were found and explained by Maslov 
and Sneppen Q on the level of interacting proteins and 
genetic regulatory interactions. According to their re- 
sults links between highly connected nodes are systemat- 
ically suppressed, while those between highly-connected 
and low-connected pairs of proteins are favored. In this 
way there is little cross-talk between different functional 
modules of the cell and protection against intentional at- 
tacks, since the failure of one module is less likely to 
spread to another one. Also in immunological networks 
one speaks of lock- and key-interactions between molecu- 
lar receptors and antigenic determinants 0. In general, 
complementarity is essential for pattern recognition inter- 
actions, underlying biological and biochemical processes 
as well as for symbiotic species in ecological networks. 
Of course, it is not at all obvious or necessary that com- 
plementarity in "internal" (functional) properties should 
be manifest in topological features like the degree-degree 
correlations. Therefore we study the relation between 
self-similarity and degree-assortativity in this paper. 



II. MEASUREMENTS AND RESULTS 

For the genetic regulatory networks Saccharomyces 
cerevisiae jj and Escherichia coli Q we observe a power- 
law behavior of N B /N for l B > 3 with d B = 5.1 ± 0.3 
for S. cerevisiae and d B = 3.4 ± 0.2, respectively. The 
obtained degree-distributions are scale-free and satisfy a 
power-law with exponent 7 ~ 2.6 for S. cerevisiae and 
7 ~ 2.1 for E.coli. The scaling relation between the 
exponents d B and dk is also satisfied within the error bars 
for both networks. 

In Table[I]we summarize the results also for some addi- 
tional networks, for which we list their properties of self- 
similarity and disassortativity. If we confirm the property 
of self-similarity it means not only the scaling behavior of 
N B /N according to a power-law, but also the numerical 
verification of the scaling relation of Eq^and the invari- 
ance of 7 under renormalizationpj. This is more conclu- 
sive, because it is sometimes difficult to disentangle expo- 
nential from power-law behavior of N B /N for networks 
with a small diameter (for example see Fig^for the reg- 
ulatory network of S. cerevisiae with an inset that shows 
the same data points on a log-linear scale instead of the 
log-log scale), whereas the scaling relation only holds for 
a power-law of the decay, it is easier to prove or disprove. 
A confirmation of (dis)assortativity refers to histograms 
with (negative) positive slope of next-neighbor degree- 
degree correlations and/or a (negative) positive Pearson 
coefficient, respectively. In most cases we measured the 
degree-degree correlations also between nodes at distance 
d = 2, 3 as indicated in the figures. 

The first four networks of Table [I] refer to the genetic 
regulatory networks of S.cerevisia and E.coli, the scien- 
tific collaboration network, and the internet on the au- 
tonomous systems level. For these networks the proper- 
ties of column 2 and 3 were examined by us, while for the 




FIG. 1: S. cerevisiae: (a) Normalized number of boxes Nb as 
a function of linear box size Ib to read off ds (b) Rescaling 
factor s(£b) as function of the box size to read off dk 



TABLE I: Networks from datasets characterized by proper- 
ties of self-similarity and disassortativity. 



Network 



genetic reg. S.cerevisia \2] yes FigHT" 

genetic reg.E.coli [£] yes Figl3l 

scient.collab. \8j no Fig[^]a 

internet aut.sys.[£] yes 
biochem. pathway E.coli [8] yes fl] 
actor [8] yes fl] 

www \8j yes [LJ 



Self-similarity Disassortativity 



yes Fig0]a 
yes Fig0]b 
no Fig b 

yes HO] 
yes [10] 
(yes, A: > 1000)Fig|U 
yes [10] 



last three networks (the biochemical pathway network of 
E.coli, the actor network and the world- wide- web) , the 
self-similarity was established before pj, and we studied 
their property of disassortativity in addition. In partic- 
ular the actor network deserves some further comments. 
The actor- network is self-similar pj , but its positive Pear- 
son coefficient suggests that it is assortative, in contrast 
to all other self-similar networks we have studied so far. 
A closer look at its next neighbor-degree-degree correla- 
tion (Fig|2J) shows an assortative behavior for degrees up 
to the order of 1000, but slowly decays for larger degrees 
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FIG. 2: Degree-degree correlation (k n (d)) for d — 1 against 
k for the actor network, showing initially assortative behavior 
for k -< 1000 
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FIG. 3: Same as Fig0 but for E.coli to read off (a) d B (b) 
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FIG. 4: Degree- degree correlation {k n (d)) versus degree k for 
distances d — 1,2,3, (a) S.cerevisiae and (b)E.coh 



and becomes disassortative. The degree-degree correla- 
tion between nodes at distance larger than 1 is decreasing 
with degree k for all k. 

Moreover, some comments are in order to the yeast- 
genetic regulatory network with 3456 nodes and 14117 
edges, (cf. Fig^). Since it has a diameter of 9, the largest 
£b value for the tiling is 10. Therefore we have only 8 
data points available for the fit. Each point corresponds 
to an average over 100 tiling configurations. Different 
tiling configurations result from different starting seeds as 
well as the random selection of neighbors during the tiling 
process. The data point at Ib = 3 in FigQ]lies clearly 
outside the fluctuations about the average over different 
tiling configurations, thus outside the error bars, which 
are at least two orders of magnitude smaller than the re- 
spective value of Nb/N, so that they are not visible on 
the scale of the figure. In Fig. [31 we find a similar behav- 
ior for E. coli. The deviation from the power-law behavior 
at Ib — 3 goes along with an assortative degree-degree 
correlation between nodes at distance d = 2 as it is seen 
from Fig^Ji and FigQJ), showing the degree-degree corre- 
lation of S.cerevisiae and E.coli, respectively, at distances 
d = 1,2,3. The data in Fig^explicitly show the disassor- 
tative behavior at d = 1 and d = 3 for both S.cerevisiae 
and E.coli. However, for d = 2, we find that there is 
a certain value of k, k = k*, at which (k n ) abruptly in- 
creases and slowly decreases for k > k* . Here k* ~ 30 
for S.cerevisiae and k* ~ 10 for E.coli. These mixed 
properties of assortativity and disassortativity seem to 
go along with the deviation from the power-law behav- 
ior of Nb/N. On a qualitative level, this is plausible if 
we focus on a hub that should be present in a scale-free 
network. In an assortative network (assortative say at 



FIG. 5: Scientific collaboration network (a) exponential de- 
cay of Nb /N(£b) (b) degree-degree correlation {k n (d)) for 
d = 1, 2, 3 against k 
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FIG. 6: Pearson coefficient for several data sets, in particular 
the PIN-data taken from [ij . 



distance d, for example d = 2), this hub is likely con- 
nected to another hub within the distance d. If this hub 
is chosen as a seed of a box in a tiling with linear box size 
£b > d, we need much less boxes to cover the many nodes 
in the neighborhood of the hub than in a disassortative 
network. 

In a network which is assortative not only for a certain 
range of k, but for all k and at distances d > 1, like the 
scientific collaboration network, Nb/N actually decays 
faster than power- like for all Ib, as it is seen from the 
exponential fit of Fig|S^. 

The scaling relation between the exponents "f,dB and 
dk assumes the scale-invariance (under renormalization) 
of the degree distribution, that is the invariance of the 
exponent 7. Similarly, it is of interest how the disassor- 
tativity transforms under renormalization (as defined in 
PJ). As we see from FigEl even networks like the scien- 
tific collaboration network (Fig|3J), which are originally 
assortative, transform to more and more disassortative 
ones under iterated renormalization. (The number of 
renormalization steps is determined by the size of the 
networks, in particular by its diameter. The final step is 
achieved when the reduced network consists of just one 
node.) Therefore the transformation behavior of disas- 
sortativity seems to be an effect of the renormalization 
procedure rather than an intrinsic self-similar property of 
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the network. Similarly, we measured the transformation 
behavior of the clustering coefficient under renormaliza- 
tion of self-similar networks. As the data |Hj show, it 
is an invariant property of scale-free networks, while it 
changes under renormalization for non-self-similar ones 
like the Barabasi- Albert one [Tlj . 

To summarize, we find numerical evidence that self- 
similar scale-free networks are preferably disassortative 
in their degree-degree correlations. For biological net- 
works this result may reflect the complementarity in in- 



teractions that is observed on various levels. 
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